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Abstract 

Second harmonic generation in negative index metamaterials is considered. Theoretical analysis 
of the corresponding model demonstrated significant difference of this phenomenon in conventional 
and negative index materials. In contrast to conceptional materials there is nonzero critical phase 
mismatch. The behavior of interacting weaves is dramatically different when phase mismatch is 
smaller or greater than critical value. 
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I. INTRODUCTION 



Experimental demonstration of the phenomenon of negative index of refraction first in 

fi n n 

the microwave [l| and latter in the optical regime [2|, |3[ has stimulated growing interest in 
nonlinear properties of negative index materials j4|. This interest is motivated by specifics 
on the interaction of electromagnetic waves with negative index materials. In combination 
with a nonlinear response of the optical material to electromagnetic radiation, this iteration 
leads to a new nonlinear optical phenomena. Study of these phenomena is of considerable 
importance both for better understanding of fundamentals of electrodynamics of negative 
index materials and their applications. One of the most fundamental property of negative 
index material is an opposite directionality of the Poynting vector, characterizing the energy 
flux, to the wave vector k. On the other hand, the negative index property can be realized 
only on particukar wavelength intervals. These two features are offering a very unusual 
type of multi-wave interactions, if frequencies of interacting waves correspond to frequency 
intervals where optical material has different signs of refractive index. Multi-wave interaction 
must satisfy a phase matching condition, which is possible only when all wave vectors are 

n 

pointed in the same direction [5j. Therefore energy fluxes of the waves with frequencies 
corresponding to a negative sign of refractive index will propagate in opposite direction to 
those with frequencies corresponding to a positive sign of refraction index. 

Such effect was suggested for the first time in , which considered the particular case of 
three waves interaction - second harmonic generation. A solution of the equation describing 
second harmonic generation in the case of exact phase matching was given in 18|]. The 
feasibility of parametric amplification using three-wave interaction for compensation losses 
in negative index materials was studied in [l^]. The dynamics of interacting wave packets 
propagating in negative index materials in the case of second harmonic generation was 
considered in [l^. It was shown that in contrast to a weak intensity of pump field, at 
high intensities a second harmonic pulse can be trapped by the pump pulse and forced to 
propagate in the same direction. 

In this paper we investigate second harmonic generation in the presence of phase- 
mismatch A. This is an important case since phase-mismatch is more relevant to realistic 
experimental conditions. Additionally, it introduces two types of spatial distribution of sec- 
ond harmonic field intensity along the sample: monotonic and periodic 1 on the coordinate. 
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Both cases are considered in this paper. We also studied second harmonic generation near a 
critical phase-mismatch value, when the material becomes transparent for the pump wave. 



II. BASIC EQUATIONS 



The system of equations describing three wave interactions (one dimensional case) in a 
- medium for the s 
following form js, I?!: 



- medium for ;he slowly varying envelope and phase approximation can be written in the 



where wave numbers fcj, j = 1, 2 are defined as follows A;| = [ujj/cf e{ujj)jj,{ujj) and kj is the 
sign of the square root of n'j = e{ujj)fi{uj) and 

P^^iui) = x\uJi; Us, -uj2)A3A; exp {tzih - k^)) 

P''\uj2) = x\^2; c^3, -uJi)AsAl exp {iz{ks - k^)) (2) 

P^^{uJs) = x\ujs; uJi, uj2)AiA2 exp {tz{ki + k2)) 

For the case of second harmonic generation Eqs. ([T]) take the following form: 



oz Vu, ot / c^ku, 



dz ^^2.dt - ' c^k2^ """P ^'^^'^ 



(3) 



where A/c = 2/c^ — k2uii A^^ is the fundamental wave with frequency cu, and a.A2u) is the 
second harmonic generated in the material. We consider the case, when the refractive 
index is negative at the fundamental frequency u and is positive at the second-harmonic 
frequency 2uj. The parameter A/c plays an important role for the spatial distribution of the 
electromagnetic field along the sample. In the next section we consider both cases A/c = 
and Afc ^ 0. 
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III. CASE OF IDEAL PHASE MATCHING Ak = 



We consider second harmonic generation for continuos waves in a medium under ideal 
phase matching conditions Ak = 0. The length of the sample we assume to be L. Using the 
symmetry properties of the susceptibility tensor with respect to permutations of u and 
2u) frequencies, the mathematical model of second harmonic generation can be formulated 



in the following way 







dz 








dz 




A^(0) = 


Al A2UL) 



A2^Al 



(4) 

(5) 
(6) 



where K = 2ttx'^{2u))/c'^ = 7rx^(u;)/c^. Let us represent the complex functions and A2U} 
in terms of amplitudes 61^2 and phases 2 



A^ = d exp (iLpi) A2UJ = 62 exp {up2). 



(7) 



Substitution of Eqs. ([7]) into (jlj) and ([5]) , and separation of real and imaginary parts lead 
to the following system of equations: 

dei 



Keie2 sin (6'), 



dz 

de2 2 ■ /Q\ 
-J— = sm {&), 

dO (A \ 



(8) 



with boundary conditions: 



ei(0) = eio, e2(i:)=0. 
Here Q and k are defined as follows 

From the first two equations an integral of motion follows: 



(9) 



(10) 



el — 62 = ml = const 



(11) 
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This integral of motion corresponds to the modified Manley-Row relation. In case of second 
harmonic generation in conventional materials the Manley-Row relation is equivalent to con- 
servation of energy (ef + e"^ = const). In our case, relation fllip corresponds to conservation 
of total flux of the energy. The second integral of motion for the system reads as: 

6^62 cos {9) = 1712 = const. (12) 

The integral of motion fll2l) is consistent with boundary conditions only if cos (6) = 0. 
Taking into account that the pump wave energy decays in z, we conclude that the phase 
difference is equal to 9 = Stt/2, therefore the system of equations (jH]) can be represented as 
follows: 

dei de2 2 rio\ 

_ = _,e,e2, ^ = -e, (13) 

The solution of f lT3|) has the following form 

ei(C) = mi/ cos(mi(/- C)) 

(14) 

62(0 = "^1 tan (mi(/- ()) 

here ( = kz and / = kL. The solutions (1141) unknown parameter mi, is the value of 
the fundamental field at the end of the sample. This parameter can be found from the 
boundary condition ([9]). Taking into account the Manley-Row relation (ITTl) . it leads to the 
transcendental equation for mi: 

eio = mi/ cos(mi/). (15) 

This equation can be solved numerically. The solution of ( lT5l) together with ( IT^ determines 
the field distribution along the sample. The dependence of intensities and on ( is 
represented in Fig{T]were the intensity boundary value e^(0) is chosen to be e^Q = 3.5, here 
/ = 1 and mi = 1. 

The solution of transcendental equation ( ITSl) for / = 1 is shown in Figj2l This plot 
illustrates the dependence of the output field intensity ei(/) = mi, as a function of eio (the 
amplitude of the fundamental field pumped into the medium) . As shown in FigJJl the formal 
solution of equation (fT5l) has multiple branches. However, only the lower branch presented 
by a solid curve has physical meaning. Upper brunches represented by dashed curves are 
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FIG. 1: The dependence of the intensity of fundamental wave ef (sohd curve) and second harmonic 
62 (dashed curve) on the distance ( with efg = 3.5 

originated from periodicity of the cos function in (fT5|) . Both ei{() and ei{() corresponding 
to these branches have singularities on the interval < C < / which is inconsistent with 
conservation of energy. Note that the lower "physical" branch shows saturation of output 
power of the electric field at the fundamental frequency ei(/) with increase of input power 
ei(0). This indicates that with the increase of input power ei(0) above 2, all excessive energy 
of pump signal converts to energy of the second harmonic signal (see Fig. |3]). 
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FIG. 2: The dependence of the intensity of output fundamental wave ei{l) on the eio 





6 



1.0 




0.8 



0.6 



ei{0)leiQ 



0.4 



0.2 



0.0 



1 2 3 4 5 6 

^10 







FIG. 3: The dependence of conversion efficiency of the pump to second harmonic fields. 

IV. SECOND HARMONIC GENERATION IN THE PRESENCE OF A PHASE 
MISMATCH A/c / 

Let us consider the impact of phase mismatch Afc (Eq. ([3])) on second harmonic genera- 
tion. The system of equations describing the spatial distribution of field amplitudes ei^2{z) 
and phase difference 6{z) in the presence of phase mismatch reads: 



Here 6 = ip2 — 2ipi — Akz and k, is defined in (|TOl) . By introducing variables ( = kz and 
/ = kL Eqs. (fT6|) can be represented in the following form: 



dei 

dz 
de2 

dz 



Keie2 sin (9), 



(16) 





6162 sin (9), 



e\ sin {9) 



(17) 




here A = I\kjK. The Manley-Row relation in this case remains unchanged: 



— 63 = = const, mi = ei(Z). 
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and a second integral of motion in presence of phase mismatch reads: 

626^ cos (9) H — ^ = 1712 = const (18) 
Taking into account boundary condition 62(1) = 0, we conclude that m2 = and therefore 

The function f|T9|) has an extremum at 63 = ml and cos (6) at this value of 62 gives cos (6*) = 
—A/Arrii. Since | cos(^)| < 1, then there is the critical value of mismatch |Acr.| = 4mi such 
that max | cos^^| = 1. Notice that flTI?]) is defined for arbitrary values of 62 if |A| < Ami. If 
|A| > Ami then there is a forbidden gap for values of 62'- 

^ (|A| - - Al) < 62 < ^ (|A| + v/A2 - Al) . (20) 
In this case | cos0| < 1 if 

e2 > ^ (|A| + VA2 - Al) (21) 

< 62 < ^ (|A| - v/A2 - Al) . (22) 

Since the value of 62 on the right side of the sample is ser to be zero (62(0 = 0), the branch 
of 62 values ( 12T]) is not accessible. Values of 62 in this case remain within the branch (l22ll . 
In this case the conversion efficiency of the pump wave to second harmonic is limited by 



the value 4eio/ (I A| — a/A^ — 16mf ). The dependence of f{9) = —Ae2/2{m\ + e^) on 62 for 
different values of mismatch is shown in Fig. HI The bold curve on this figure corresponds 
to a critical value of the mismatch. The forbidden gap for 62 can be seen for two lowest 
curves, when curves are below —1. The presence of a forbidden gap for 62 suggests the 
existence of two types of solutions for 62. The first type corresponds to mismatch values 
|A| < Ami and in this case 62 is not bounded from above. This means that the conversion 
rate of fundamental harmonic to second harmonic, in principle, can be high (close to 1 - 
ideal conversion, similar to Fig. |3] in the previous section). The second type of solutions 
correspond to mismatch values |A| > 4mi; in this case the amplitude 62 is bounded from 
above < 4e2 < |A| - v/A^^OGmf. This means the is a limitation of the output intensity 
of second harmonic field at the growing input intensity of fundamental harmonic. 

For further considerations it is more convenient to deal with field intensities rather then 
with amplitudes. Using expression ( 1T9|) for cos (6'), the second equation of (fT7|) can be 
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FIG. 4: The dependence of the cos [9) on 62 with different values of A = /c x mi, here k = 1, 2, . . . , 6 
and mi = 1.2. The bold curve corresponds to a critical value of A = Ac^ 



represented as an equation for the intensity P2 = e\. 



dc 



1/2 



(23) 



where F[P2) is a quartic polynomial 

F(P2) = 4P| + (8m2 - A^) P| + 4m^P2 

with the following roots: 

1 



P2c = ^ ( - 8m^ + A^/A2 - 16m2 
1 



P26 = ^ ( A^ - 8m^ - A./A2 - 16m? 



P2a = 0. 



(24) 



Notice that these roots ( |2^ define the forbidden gap [a/^, a/Pjc] for values of 62 (see Fig. H] 
and equations (12T]). (I22l)). 

Based on this qualitative analysis, we conclude that there are three regimes of second 
harmonic generation controlled by the absolute value of the phase mismatch. In the following 
subsection we will analyse solutions describing spatial field distribution inside the sample. 



A. Three regimes of second harmonic generation 

The absolute value of the phase mismatch determine three different regimes of second 
harmonic generation: |A| < A^r, |A| = A^^ and |A| > A^^. First we consider the case of 
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FIG. 5: The dependence of the output field ampUtude ei(Z) at the fundamental frequency on eio 
near the critical value of phase mismatch Acr 

subcritical mismatch: |A| < Acr- 

1. Subcritical mismatch 

In the case where |A| < Acr, roots are complex- valued and the solution of fl2^ can 



be expressed in terms of Weierstrass function p 2l|. By expanding F (P2) into Taylor series 
and introducing a new variable: 



4(P2 - P2a) 24 

Here derivatives of the polynomial are taken with respect to P2, then the solution of equa- 
tion fl23l) can be represented in an implicit form: 

f°° ds 
Js [As-^ - g2S - gs} 

where g2 and g^ are invariants of Weierstrass function. 

92 = ^ {8ml -A'Y- 4m{, g, = ^m^ (8m? - A^) - ^ (8m? -A')'. 
Finally, the amplitudes of second and fundamental harmonics have the following form: 
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ei(C) = ^mf + eKC), (26) 

e2(C) = = (27) 

V(K^-C;fe^73)-(8m?-A2)/12) 

The parameters g2 and (73 are functions of A and nii. To determine solutions of Eqs. f ll7p 

we need to solve for the unknown value of the output pump wave mi. The value of mi can be 

found taking into account the output boundary condition mi = ei(Z) and the Manley-Row 

relation (fTT]) . which lead to the following transcendental equation for mi: 

4 

= + p(/;^2,^3)-r8m?-A^)/12 (2^) 
To determine the unknown parameter mi = ei(/), Eq. (!28!) needs to be solved numerically. 
The analysis of the ei(Z) dependence on eio and A shows that with increasing phase mis- 
match from to Acr, all branches (physical and nonphysical) shift upwards and nonphysical 
branches change their shapes. Fig. [5] shows the dependence of the output field amplitude 
ei(/) at fundamental frequency on eio near the critical value of phase mismatch A^r- Sheet, 
labeled as "1", corresponds to the physical branch, while sheet labeled as "2" represents the 
first nonphysical branch. Other nonphysical sheets are located above nonphysical sheet "2" 
shown on Fig. O The spatial distribution of ei(C) and 62 (C) can be found by substitution of 
the solution of the Eq. ( !28l) (mi) in Eqs. ( l26l) and ( l27l) . We found that all solutions ei(C) 
and 62 (C) are monotonically decreasing in C,. An example of 62 (C) at A = Acr/2 is shown in 
Fig. ED. 

The conversion efficiency a = e2(/)/eio as a function of the input amplitude eio is pre- 
sented in Fig. [71 As one can observe, a is approaching its asymptotic value a = 1 in slower 
fashion for larger values of |A|. 



2. Critical mismatch 

When the value of phase mismatch is critical |A| = A^, the roots P2C = P2b = ^Irit/^^ = 
ml and the discriminant of the Weierstrass function is zero. In this case, the function 
p{l — (192,93) can be represented in terms of hyperbohc functions. Thus the Eqs. fl26|) 
and (127]) take the form: 

ei(C) = ^Jmj + el{0 (29) 
62(0 = mitanh(mi(/- C)) (30) 
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FIG. 6: The dependence of second harmonic's amphtude 62 (C) on the Q with different values 
of phase mismatch(the solid curve: A = 0, large dashed curve A = 2?tii, small dashed curve: 
Acr = 4mi, intermediate dashed curve: A = lOmi , doted curve: A = 25mi). 

and the transcendental equation for mi reads as 

e^o = ml{l + tanh^ (mi/)) (31) 

The numerical solution of Eq. ( |3T]) is shown in Fig|5] (line "3" ) . Observe that at large values 
of eio, the solution of Eq. ( l3Tl) is proportional to eio {rrii ^ eio). Therefore, at large values of 
eio the conversion efficiency a — )■ tanh(Zeio) is always less then one while in the subcritical 
regime a — )■ 1 (see Fig. [7]). 




3. Overcritical mismatch 

At large mismatch values, when |A| > A^., roots of are real. In this case it is 
convenient to represent p{l — (',92,93) in terms of Jacobi elliptic sn function {2^. Eqs f l26|) 
and (1271) can be represented as 



e2(C) 



P2h sn 



here 7 = ^/P2b/P2c, and the equation for mi takes the following form: 



^10 = "^1 + P2b sn^ 



'P2cl,l 



(32) 
(33) 

(34) 



The sheet corresponding to solutions of (1M|1 is labeled in Fig. l5l as "4" . In contrast to the 
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FIG. 7: Dependence of conversion efficiency a = e2(0)/eio on input field amplitude eio with 
different values of phase mismatch A: Solid curve A = 0; large dashed curve A = 3.5mi; dashed 
curve A = 3.9mi; small dashed curve Agj. — Attix] doted oscillation curve A — 4.5?7ii > A^r 

subcritical regime, all solutions in this case are represented by a single sheet. This sheet 
has folds. Hence the intersection of this sheet with plane corresponding to A = const gives 
multivalued dependance of ei(/) on eio- This dependance for two different values of A is 
shown in Fig. [8l 

In the supercritical regime the second harmonic field experiences spatial periodic oscil- 
lations with period 4K(7) (see Fig. [6]). The distance between neighboring zeros C of the 
amplitude of second harmonic is determined by the following formula: 



If the slab length satisfies the condition / = n x C, {n = 1,2,3...), then the amplitude of 
the second harmonic wave is zero at the both ends of the slab (zero conversion efficiency). 
Therefore such slab is transparent for a pump wave. A plot of the transmission coefficient 
53 = ei(/)^/efQ as function of eio is shown in Figj9l The transmission coefficient is equal to 
1 at the points labeled as "1", "2", . . . (transmission resonances). The spatial distribution 
of the fundamental and second harmonic fields corresponding to the transmission resonance 
at the point "1" (see Fig|9]) is shown in FigJTOl 



c 



2K(7) 



(35) 
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FIG. 8: The dependence of the output field amphtude ei(Z) at fundamental frequency on eio- Solid 
curve A = 4.1mi, dashed curve A = 4.5mi 

V. CONCLUSION 

We considered second harmonic generation in negative index materials. Specifics of this 
process is in the negative value of refractive index for the pump wave and the positive value 
for second harmonic. This led to important features which are different from the case of sec- 
ond harmonic generation in conventional dielectrics. The main difference is in the existence 
of nonzero critical values of the phase mismatch. If the absolute value of phase mismatch is 
below critical, then the field intensities are monotonically decaying along the sample. When 
the absolute value of phase mismatch exceeds a critical value, monotonic decay of intensities 
transforms to a spatial periodic oscillations. Note, that in the conventional case the critical 
value of phase mismatch is zero. 

Another important feature is the dependance of conversion efficiency on the amplitude 
of the incident pump wave. When the absolute value of phase mismatch is below critical 
value, then the conversion efficiency asymptotically approaches 100% at large values of the 
incident pump wave amplitude. It should be stressed that in this case the asymptotic value 
of conversion efficiency does not depend on the phase mismatch value. The phase mismatch 
affects only the rate of approaching of conversion efficiency to its asymptotic value. When 




14 




2 4 6 8 10 12 14 

ei(0) 

FIG. 9: Dependence of the transmission coefficient S> on pumped field amplitude eio- A = 4.2mi 
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FIG. 10: Spatial distribution of the intensities ef{Q (solid curve) and 63 (C) (dashed curve) inside 
the slab 

the phase mismatch is exactly equal to critical value, then the asymptotic value of conversion 
efficiency experiences a jump to a value which is less than 100%. When the absolute value 
of phase mismatch is above critical value, the conversion efficiency becomes an oscillatory 
function of the incident pump wave amplitude. 

Finally, we found that the dependance of output amplitude of the pump wave on its 
input amplitude is single valued if the absolute value of phase mismatch is below critical 
and becomes multi- valued in the opposite case. 
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